Laser-driven nanoplasmas in doped helium droplets: 
Local ignition and anisotropic growth 
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Doping a helium nanodroplet with a tiny xenon cluster of a few atoms only, sparks complete ioniza- 
tion of the droplet at laser intensities below the ionization threshold of helium atoms. As a result, 
the intrinsically inert and transparent droplet turns into a fast and strong absorber of infrared light. 
Microscopic calculations reveal a two-step mechanism to be responsible for the dramatic change: 
Avalanche- like ionization of the helium atoms on a femtosecond time scale, driven by field ionization 
due to the quickly charged xenon core is followed by resonant absorption enabled by an unusual 
cigar-shaped nanoplasma within the droplet. 

PACS numbers: 87.15.ht,31.70.Hq,36.40.Gk,36.40.Wa 



Helium nanodroplets are fascinating objects due to their 
multi- facet properties, ranging from superfluidity itself 
[li] to unusual electron dynamics upon photon impact [3] ■ 
Most widely, however, helium droplets are used as a "cat- 
alyst" to facilitate an intended physical process without 
altering it. A prime example is the cooling of molecules 
inside a droplet for high-resolution spectroscopy [1, 0], 
the creation of unusual molecules on the droplet surface 
, but also the assembly and transport of clusters inside 
a helium droplet . The catalytic property originates in 
the high ionization potential of helium in connection with 
the almost frictionless, superfluid environment. 

Here, we will show that also the opposite is possible: 
By doping the droplet with a handful of heavier rare- 
gas atoms the inert and transparent helium droplet can 
be turned into a highly reactive object which absorbs 
infrared light very effectively. In this case, the heav- 
ier atoms serve as a "catalyst" to activate the helium 
droplet. This finding is quite surprising since the pris- 
tine droplet cannot be ionized at all with laser light of 
780 nm wavelength at an intensity of / '--^ 10^'*W/cm^ 
which we will apply. Yet, with a few xenon atoms in- 
side, all electrons from the helium atoms are removed so 
that the entire droplet containing as many as 10^ helium 
atoms turns into a nanoplasma. The extremely efficient 
ionization - even for large droplets - requires two ele- 
ments, an initial seed and a resonant energy absorption 
process from the laser pulse. The latter occurs on an elec- 
tronic time scale, in contrast to the well known resonant 
absorption during the Coulomb explosion of a homoge- 
nous cluster which occurs on the much slower time scale 
of nuclear motion 0, Q . 

In order to take the role of a seed the ionization 
potential of the embedded species must be lower than 
that of helium (i?ip=24.6 eV). We have used xenon 
(£'ip=:12.1 eV) for the present investigation. As known 
from quantum Monte Carlo calculations 0, rare-gas 
dopants settle in the droplet center. When the laser pulse 
is ramped up, xenon is ionized first producing a strong 
static electric field, which is capable to significantly ionize 



helium in combination with the laser field. In this first 
phase, electrons from the droplet migrate to the center 
of the cluster, but also leave the cluster, i.e., are outer 
ionized, increasing the total cluster charge and thereby 
enhancing helium ionization. However, the degree of ion- 
ization per helium atom, which can be achieved with field 
ionization, decreases rapidly with the helium droplet size, 
since the static field of the ionic charge falls off quadrat- 
ically with the distance from the center. 

Hence, to turn very large droplets into completely 
stripped atoms, another, even more powerful mechanism 
for ionization has to kick in and this is a new kind of 
plasma resonance which leads to an ionization avalanche. 
A little thought reveals that resonant plasma absorption 
should not be possible because the dipole plasmon fre- 
quency of a spherical helium nanoplasma is too high to 
come into resonance with a short laser pulse of frequency 
Lo. However, since the ionization starts from a small seed 
at the droplet center and is driven by a linearly polarized 
laser, a non-spherical, cigar-shaped nanoplasma forms. It 
has a plasmon frequency along the longer axis lower than 
that of a sphere with an equal charge density. Therefore, 
resonance with the laser frequency is possible and hap- 
pens in fact very fast, within a few femtoseconds. As a re- 
sult, helium droplets containing more than 10^ atoms are 
completely ionized with the help of 13 or fewer embedded 
xenon atoms, while in the absence of the xenon seed, the 
droplet is not ionized at all. The ionization avalanche can 
be invoked by a wide range of laser frequencies, intensi- 
ties and pulse durations. It is also robust with respect to 
the helium density and the species of seed atoms as long 
as they have lower ionization potential than helium. 

In order to demonstrate the two elements necessary 
for the avalanche ionization and in particular the new 
anisotropic nanoplasma resonance, we will present cal- 
culations for A — 780 nm and for A = 200 nm laser 
light, because in the latter case the laser frequency is 
too high to be matched by an eigenfrequency of the 
anisotropic nanoplasma and only the first mechanism is 
active, namely field ionization of the helium droplet. We 
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FIG. 1: Ionization and energy absorption of a He25oo droplet 
doped with a Xen cluster, (a) Charge per helium atom and 
(b) absorbed energy as a function of the number of xenon 
atoms n. The charge q refers to the average number of elec- 
trons released from a droplet atom. Both quantities are shown 
for two different laser wavelengths: A — 780 nm (circles) and 
A = 200 nm (squares). The other laser parameters were iden- 
tical: peak intensity 7 = 7 x 10^'*W/cm^, Gaussian pulse 
exp(— In 2(f/T)^) with a duration T = 20 fs, linear polariza- 
tion. 



use a classical molecular dynamics approach described 



elsevifhere [10|, which we have developed previously [8l.lll| 
and which has been applied in a similar way also by other 
authors 13, 13, 14 1. For droplets with more than 2500 
atoms we have used an implementation fl5i| of the fast- 
multipole method 16] to calculate the Coulomb interac- 
tion of electrons and ions. 

Firstly, we will show that the avalanche ionization ef- 
fect occurs over a wide range of sizes of embedded xenon 
cluster n - starting from a very few atoms - and also over 
a wide range of helium droplet sizes m, spanning more 
than two orders of magnitude (Figs.[T]and[2]). Secondly, 
we will analyze the ionization dynamics of a particular 
cluster to understand the peculiarities of the resonance 
mechanism here at work. 

Figure [T] summarizes the dependence of the avalanche 
ionization on the xenon core size n for a fixed droplet con- 
taining m — 2500 helium atoms. As can be seen in Fig. [2 
already very few embedded xenon atoms lead to complete 
inner ionization, whereby inner ionization refers to the 
(average) number q of electrons released from a droplet 
atom, with the maximum value q — 2 for helium. Of 
course, these electrons may still be trapped by the global 
droplet potential forming a plasma as discussed below. 
Outer ionization, i. e., the removal of electrons from the 
entire cluster, is determined by the plasma dynamics and 
may depend on post-pulse effects 14 1. 



The average charge of two indicates that the entire 
droplet has turned into a nanoplasma, i. e. there are nei- 
ther neutral atoms nor singly charged helium ions left. 
How dramatic this effect is, becomes also evident in the 
absorbed energy, cf. Fig.[T]3. Whereas the absorption for 
one xenon atom is negligible, it reaches for cores with 
n — 4... 13 xenon atoms a few hundred keV due to 
the ignition of the droplet. The increase occurs more 
abruptly for the longer wavelength (A ~ 780 nm, circles 
in Fig.[T]) than the shorter one (A — 200 nm, squares in 
Fig.[T]) as can be seen for the charge q in panel (a) as well 
as for the absorbed energy in panel (b). The energy in- 
creases further for cores with n = 55 ... 135 xenon atoms 
up to values of about 1 MeV, whereby now the core itself 
contributes considerably to the absorption. 

The difference in absorption between the two laser 
wavelengths becomes much clearer from the dependence 
on the helium droplet size m in Fig. [21 where the size of 
the core Xei3 was kept constant. Increasing the droplet 
size m by almost two orders of magnitude does not 
change the behavior for the larger wavelength (circles 
in Fig. [2]): Up to a size of m ~ lO*^ we find complete 
inner ionization. We recall that not a single of these 
helium atoms would be ionized without the xenon core. 
In contrast, at shorter wavelength the charge per atom 
decreases for larger droplets (squares in Fig. [2]) since the 
static field generated by the xenon core ions can only ion- 
ize a droplet of a certain size as the field strength drops 
quadratically with the distance. Obviously, the ratio of 
xenon to helium atoms is crucial. For smaller droplets 
m ^ 10^ the core ions alone drive the complete helium 
ionization, for larger droplets m » 10'' they do not. The 
reason for the different behavior at the two frequencies 
is resonant absorption. We will discuss it in detail since 
it exhibits novel features which do not occur in homo- 
nuclear clusters. 

Helium's nuclear charge of two does not allow for 
charge densities q larger than two times the particle den- 
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FIG. 2: Charge per helium atom as in Fig.[l^ for helium 
droplets doped with a Xeia core as a function of the droplet 
size m. The laser parameters are the same as in Fig.[T] The 
black circle marks the droplet shown in Fig. [31 
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FIG. 3: (color online). Anisotropic nanoplasma forma- 
tion, (a) Extensions Ly and L± of the nanoplasma parallel 
(blue/upper line) and perpendicular (red/lower) to the laser 
polarization as a function of time t for the droplet marked 
in Fig. [2] (m = 20000). Lines are averaged over a laser period 
in order to smooth sub-cycle structures. The gray-shaded 
area sketches the laser pulse envelope, (b) Ratio Ly/Lx as 
obtained for the smoothed curves from panel a. 



sity of the droplet: ftnax = 0.04/A'^. The correspond- 
ing eigenfrequency for a dipole oscillation of an electron 
cloud with respect to the ionic background, that has a 
uniform charge density g, reads for spherical geometry 
" with the plasma frequency flp\ = \f^nQ [17j, 



Sect. 2.1.2] resulting in Omax — 0.16 a. u. for a completely 
inner- ionized helium droplet. With this eigenfrequency it 
is impossible to match the laser frequency lo = 0.23 a. u. 
for the shorter wavelength A = 200 nm. On the other 
hand, ilmax is considerably larger than uj = 0.058 a. u., 
the laser frequency for A = 780 nm. To reach the reso- 
nance, the eigenfrequency of the droplet could be lowered 
by expansion of the sphere, which is the case for the usual 
nanoplasma resonance 0, This, however, can be ruled 
out here, since we observe resonant absorption for a pulse 
length of 20 fs which is too short for sufficient droplet ex- 
pansion, even for the light helium nuclei. In fact, the 
nuclei are almost frozen during the pulse. Therefore, the 
eigenfrequency of the droplet must be lowered by a dif- 
ferent mechanism to come into resonance with the laser 
frequency. 

The key to understand this mechanism is the geome- 
try of the nanoplasma formed. Our calculations reveal 
that ionization starts from the doped core in the cen- 
ter of the droplet and subsequently spreads towards the 
droplet surface. Most importantly, this spreading oc- 
curs anisotropically, namely faster along the laser po- 
larization than perpendicular to it, as shown in Fig.[3l 
There, we have plotted the 3-dimensional extension of 
the nanoplasma using 



and corresponding definitions for Ly and , where qi {t) 
is the charge state of the ith atom located at — 
(xi,yi, Zi). The sum runs over all atoms. Since our sys- 
tem is axially symmetric for a laser linearly polarized in 
z-direction, we plot the extension iy|| = along the laser 
polarization z and perpendicular to it, Li_ = [Lx+Ly] /2. 
Figure [3^ shows clearly that the nanoplasma spreads 
much faster along the laser polarization axis with lit- 
tle extension perpendicularly to it until the peak of the 
laser pulse at t = is reached. After the laser peak 
for times t > the plasma grows with increased rate 
perpendicularly to the polarization. At the same time, 
the extension along the polarization axis decreases since 
more atoms close to the "equatorial plane" with z « 
are ionized. From about 15 fs on the whole droplet is 
ionized giving rise to a spherical nanoplasma for which 
— ^||- this time all helium ions are charged and 
the extensions characterize the droplet size. The slight 
increase of both values for t > 15 fs reflects the onset of 
droplet expansion. In short, the ionization starts from 
the few xenon atoms in the center of the droplet form- 
ing an axially symmetric plasma channel. This plasma 
channel gradually elongates, cf. upper inset of Fig.lH and 
reaches the surface of the droplet. Only after that, the 
ionization spreads into the whole droplet. 

To get an estimate on the plasma frequencies involved, 
we approximate the elongated plasma channel by a cigar- 
shaped axially symmetric ellipsoid with an extension 
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FIG. 4: (color online). Elongated nanoplasma geometry. 
Eigenfrequency for an ellipsoidal nanoplasma in the di- 
rection of the long axis z in terms of the eigenfrequency for 
a spherical nanoplasma f2 as a function of its aspect ratio 
a = R\\/R± according to Eq. ([Sjl. Upper inset: Cross sec- 
tion of the nanoplasma (red/gray atoms) around the xenon 
core (blue/darkgray atoms) in the helium droplet (lightgray 
atoms). The arrow marks the laser polarization. Lower 
insets: Cross sections of equipotential surfaces in homoge- 
neously charged ellipsoids with aspect ratios of a=2 and q=4, 
respectively. 
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= Rc?l^ along the laser polarization and a radius 
i?^ = Rja^l'^ perpendicular to it. Hereby, R is the ra- 
dius of a corresponding spherical nanoplasma with the 
same volume and a = R\\/R±_ > 1 the aspect ratio of 
the ellipsoid axes. The potential inside a homogeneously 
char ged ellipsoid follows from the general integral expres- 



sion 



18j, § 99] and reads 



cl){p, z; a) = 7TQ[[l~g{a)] + 2g{a)z^ - h{a)R^\ , (2) 

with g{a) = [gin (a+/(a)) //(a) - 1]/ P{a), h{a) = 
[l+P{a)g{a)\ /a^/s, and /(a) = Vc?±X_ the limit of 
equal axes R\\ = R± the potential Eq. ^ reduces to that 
of a sphere (j){p, z; a^l) = 2ttq [{p^+z'^) /3 - R^] . Note 
that the potential Eq. ([2]) is harmonic in p and z with 
ellipsoids as equipotential surfaces shown in the lower in- 
sets of Fig. 01 However, the shape of these ellipses does 
not coincide with the plasma shape, and therefore the 
ratio of the eigenfrequencies differs from the aspect ra- 
tio of the nanoplasma. The separation of p and z in 
Eq. ^ makes it straight-forward to calculate the eigen- 
frequency of a collective dipole oscillation along the z 
axis, i. e. along the laser polarization. We get for this 
frequency ^f{a) = ^(j>{p, z;a)\^^_^^^ from Eq. Q in 
terms of the spherical eigenfrequency defined above 



(3) 



which is shown in Fig.[H For all values a > 1 the eigen- 
frequency is lower than Q, with a reduction of almost 
50% for an aspect ratio of a = 3. 

This explains the resonant energy absorption: As the 
plasma channel elongates due to the inner ionization 
along the polarization direction, the aspect ratio grows 
and the eigenfrequency decreases. At aspect ratios of 
a « 3 . . . 4, cf. Fig.[3|3, the ellipsoidal eigenfrequency is 
two times smaller than the corresponding eigenfrequency 
n of a sphere. Around this aspect ratio, f^n coincides with 
the laser frequency giving rise to the plasma resonance. 
As a consequence, enough laser energy is absorbed by the 
plasma electrons, that the ionization continues in the per- 
pendicular direction and encompasses the whole droplet. 
This happens very fast, since the resonance condition in 
form of the suitable aspect ratio of the nanoplasma is 
created by an electronic process. It does not rely on the 
slower nuclear motion of an exploding droplet as it is the 
case in the conventional resonant heating of a cluster. 

To summarize, we have shown that helium nano- 
droplets doped with small xenon clusters turn into 
nanometer-sized plasmas with all helium electrons re- 
moved. This can be achieved with a laser pulse of mod- 
erate intensity for which the pristine helium droplet is 
fully transparent. We observe an ionization avalanche on 
a fast timescale of a few femtoseconds. The phenomenon 
is robust over a wide range of laser pulse parameters and 
could be experimentally seen by detecting the depletion 



of droplets from the beam or by measuring the abun- 
dance of helium ions provided the droplets were doped 
by small xenon clusters. 

While we have discussed it from the perspective of the 
helium nanodroplet, one can also view it in a more gen- 
eral context of composite clusters, where interesting elec- 
tron dynamics has been observed for so called core shell 
systems (Ar-Xe clusters) under ultraviolet laser pulses 
[l9|. In another experiment, a strong enhancement of X- 
rays from laser-irradiated argon clusters doped by a few 
percent of water molecules was observed . We have in- 
dications that this can be explained by a seed cluster sim- 
ilarly as in the present context of the helium droplet 21 1 . 
In another scenario, deuterons, fast enough to induce nu- 
clear fusion 22] can be generated in hetero-nuclear clus- 
ters [H, certainly also spectacular as the complete 
stripping of helium in large nanodroplets sparked by a 
few xenon atoms as discussed here. 

Finally, the ionization avalanche effect in non-spherical 
plasmas which we have discovered may be of interest 
far beyond cluster physics. Similar phenomena should 
also occur for laser-illuminated solids doped with easy- 
ionizable impurities. This could be a way to create and 
study microscopic plasmas inside solids [24 1 . 
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